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CN ■ Abstract 



This paper addresses an open problem recently posed by V. Kozlov: a rigorous 
proof of the non-integrability of the geodesic flow on the cubic surface xyz = 
1. We prove this is the case using the Morales- Ramis theorem and Kovacic 
algorithm. We also consider some consequences and extensions of this result. 



1. Introduction 



In two recent papers [6^ and [f\ Kozlov posed the following open problem: to 

rigorously prove the non-integrability (in the sense of Louiville) of the geodesic 

^ ' flow on the surface xyz ^ 1. In what follows we will exploit the Hamiltonian 

nature of the geodesic equations by examining the variational equations about 
.^b ■ a planar geodesic. The crucial theorem we shall make use of is due to Morales- 

^SJ ! Ruiz and Ramis, which we quote from Q: 

\- j ' Theorem (Morales- Ramis) . For a 2n dimensional Hamiltonian system as- 

p^ , sume there are n first integrals which are meromorphic, in involution and inde- 

pendent in the neighbourhood of some non-constant solution. Then the identity 
component of the differential Galois group of the normal variational equation 
. ^^ (NVE) is an abelian subgroup of the symplectic group. 

^^ 
; I I On two-dimensional manifolds the normal variational equation (which we 

Cu ' shall derive in the next section) is a second order linear ordinary differen- 

tial equation with meromorphic coefficients. In particular, we will see that 
the NVE's of interest are Fuchsian. To show that the geodesic flow is not 
meromorphic-integrable it suffices to show that the NVE is not solvable in the 
sense of differential Galois theory: the identity component of the differential Ga- 
lois group of the NVE is not abelian. This means we cannot "build" the solutions 
from the field of meromorphic functions by adjoining integrals, exponentiation 
of integrals, or algebraic functions of elements of the field of meromorphic func- 
tions. To test this we make use of the Kovacic algorithm. 
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Before we state the Kovacic algorithm we note that this algorithni is very 
robust and can treat any second order hnear ODE with rational coefficients, 
however in the present work we will only need a very limited portion of the 
algorithm. Thus to save space we will present a much abbreviated version and 
refer the reader to the original article of Kovacic [5| and the reduced version 
appropriate for Fuchsian ODE's in Churchill and Rod [J], whose notation we 
will follow most closely below. 

Consider a linear ode of the following form 






■j^=e=r{y)t r{y)€C{y). (1) 



If the equation is Fuchsian, that is it admits only regular singular points, then 
we can decompose r{y) as 



3 



where k is the number of finite regular singular points at locations y = aj. When 
X] iJj = then y = oo is also a regular singular point, with /3oo = J2Wj + ^j'^i)- 
The indicial exponents are 



T± = i l±yTT4^ , r± = i l±Vl + 4/3, 



at 2/ = Uj and y = oo respectively. 

Kovacic proved in [5| that there are only 4 possible cases for the differential 
Galois group of ([T]) . We will see in Section 3 that we can rule out two of these 
cases immediately (cases / and ///), and as such we will present only necessary 
conditions for these cases. 

Theorem (Kovacic). Let Q he the differential Galois group associated with 
([1]), and note Q C 5*^(2,0). Then only one of four cases can hold: 

(I) Q is triangulisable (or reducible) . A necessary condition for this case to 
hold is that, defining the 'modified exponents ' a as 

1 if Pj = and Sj ^ 0; af ^ if Pj = 5^ = 0, 
= z//3oo = 0. 
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(II) Q is conjugate to a subgroup of the 'DP' group, in the terminology of 
Churchill and Rod. A necessary and sufficient condition for this case to 



hold is that, defining the following sets 

Ej = {2 + e^/TTWj, e = 0, ±2} DZif (3, ^ 0; 

E, - {4} if p, = 0, 5, ^ 0; E^ = {0} if P, = S, = 0; 

Soo = {2 + eVl + 4/?oo,e-0,±2}nZ2//3oo^0; 

^oo = {0,2,4} i//3oo=0, (2) 

there is some combination of ej e Ej and Coo € Eoo, not all even integers, 
so that 

k 

j=i 

and there exists a monic polynomial P{y) of degree d which solves the 
following ODE: 



P'" + idP" + (36'2 + W - ^r)P' + (61" + 36*61' + 6'^ - 4r6' - 2r')P = 0, 



2 ^ y-c 



(3) 



(77/j tj is finite. A necessary condition for this case to hold is that all indicial 
exponents t- and r^ are rational. 

(IV) Q = SL(2,<C), whose identity component is not abelian and therefore ^ 
is not solvable. 

In the next section we will derive the NVE about a planar geodesic on the 
Monge patch with a plane of symmetry; in Section 3 we will prove, using the 
theorems presented in this Section, that the NVE is not solvable and therefore 
the geodesic flow is not Louiville integrable on the surface xyz = 1. In Section 4 
we will consider some extensions and consequences of this result, and in Section 
5 we finish with some conclusions. 

We note that the approach followed in this paper has been used to prove the 
non-integrability of a number of problems in mechanics and celestial mechanics 
(see, for example, ^, [HI, Q, ^, fl, ^, [i|,ii[,to name but a few), but 
with the exception of another work by the author [l5| this approach is novel in 
examining geodesic flow. 

2. Derivation of the NVE 

The key feature of the surface xyz = c, where w.l.o.g. we can set c = 1, 
which facilitates this analysis is that by a simple rearrangement and rotation of 
7r/4 about the z-axis we can write z = l/(a;^ ~ y^), or more generally 

z^fi^.y), fAo.y) = 0- (4) 



This means the surface is a Monge patch (or graph) with a plane of symmetry 
(or invariant plane) , the y-z plane. The surface is actually made up of 4 identical 
components, and to demonstrate the non-integrability of the geodesies of the 
surface we need only demonstrate this property on one component. We will 
restrict our attention to the quadrant 

{a;,y,z€R3 . |^| > |y|,a;>0} 

which immediately rules out any possible divergences in /. 

To keep the approach of this section general and to facilitate the analysis of 
Section 4 we will derive the NVE on the Monge patch with a plane of symmetry 
as in dH). 

Lemma 1. The normal variational equation about the planar geodesic on the 
Monge patch Q is 

where /C|o is the Gauss curvature evaluated along the planar geodesic. 

Proof. Using the standard parameterisation {x,y, f{x,y)) and resulting line el- 
ement we may calculate the Christoffel symbols and geodesic equations: 



J^Kc^'i'^O, x'-^ix.y) 



b,c 

where a dot denotes differentiation w.r.t. arc-length s. The x — plane is 
invariant since Tyy — when f,x{0, y) — 0, and thus there is a planar geodesic 
(0, y(s), /(O, y{s))) where y solves (1-1-/^^(0, yY)y'^ = 1. Linearizing the geodesic 
equations about this planar geodesic the normal variational equation will simply 
be the variation in the a;-direction, namely 



^+{K...y 
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The problem with this equation is that the coefficient is a function of y{s), 
which is defined implicitly as a solution of (1 + /^)y'^ = 1. Clearly y will also 
parameterise the planar geodesic and thus we make the change of independent 
variable to y, calculating derivatives such as (dropping the tildes) 

d _ . d _ 1 d 



^* ' '^y Ji + n'^y 



and so on, to arrive at the NVE (dropping the subscript) 



j-.yj-.yy \ c' i ( J-yyJ-^ ^ 



^" _ ( LvlfL ) ^' + ( Jjill^ ) ^ ^ 0. (5) 



If, for a given surface z = f{x, y) with /,a;(0, y) = 0, this NVE is not solvable 
in the sense of differential Galois theory as described in Section 1, then the 
geodesic flow on that surface is not integrable. This is precisely what we will 
show in the next Section for the surface xyz = 1. Before we do however, we can 
make a comment about ([S]) : 

Notice that the equation is of the form ^" — f.yQ{y)C + f,xxQ{y)£, ~ 0. If 
f,xx/f,y — 1/y, then ([5]) would have the simple solution ^i = ciy from which we 
could construct a second solution via integrals and exponentiation of integrals 
of meromorphic functions. This leads us to consider solutions of the PDE 

yf,xx - f,y = (6) 

(where we evaluate the derivatives of / at x = 0) as candidates for surfaces 
with integrable geodesic flow. Examples include well-known integrable sur- 
faces such as f{x,y) = f{x'^ + y^) and more interesting surfaces such as / = 
cos(a;a;)e~" ^ /^. But we should not divert too much attention to (jH]): a sur- 
face which solves this equation need not have integrable geodesic flow, it merely 
passes this integrability test. 

3. Non-integrability of the surface xyz — 1 

In the case of z = f(x,y) = {x^ — j/^)^^, the NVE takes the form 



18(2 + %^), 






where we have removed the ^' term from ([SJ via the standard transformation 
[5[, and we extend the independent variable to the complex domain. There are 
8 regular singular points, aj = {0,pi, . . . ,p6, } and oo where pi denotes the 6 
roots of 2/^ -f- 4 = symmetrically distributed about the circle of radius ^a/4 
centred on the origin. We find the /3 coefficients are 

9 5 5 



Pj-[ 4' 16'- ■•'16/' ^°°-^' 



and only 5oo =^- We can see immediately that Tq = ^(1 ± i-\/8), and therefore 
case /// of the Kovacic algorithm can be ruled out (the finite case). What's 
more, none of the other r- are complex so case / of the algorithm can also be 
ruled out (the triangulisable case). 

Case // is more problematic. The sets described in ^ are 

Eq - {2}, Ei...Q = {2, 5,-1}, Eoo = {0, 2, 4}. 

There are 21 combination of the elements of these sets leading to each of d = 
and d = 1, and 1 leading to each of d = 2,3,4 (for example, d — ^(4 — (2 — 
1 — 1 — 1 — 1 — 1 — 1)) =4). For each of these combinations we must attempt 
to construct a monic polynomial of order d that satisfies ([3|) . This can be done 



using a computer algebra system such as Mathematica; the calculations are 
tedious rather than difficult. By checking each combination we can see that 
there is no polynomial P satisfying ^. We can now state the main result of 
this paper: 

Theorem 1. The geodesic flow on the surface xyz — 1 is not integrable in the 
sense of Louiville with meromorphic first integrals. 

Proof. The differential Galois group of the normal variational equation Q does 
not fall into case /, // or /// of Kovacic's algorithm, as we have shown above. 
Therefore we must have Q — SL{2,C), the identity component of which (also 
SL{2, C)) is not abelian. By the Morales- Ramis theorem of Section 1 this means 
the geodesic equations are not Louiville integrable with meromorphic first inte- 
grals. D 

4. Extensions and limitations 

It seems natural to ask can we use the same techniques to examine other 
surfaces similarly defined. We will consider two generalisations, a;"y"z" = 1 
and x"y"z = 1. 

4-1. Surfaces of the form a;"y"z" = 1 

While it might seem "obvious" that xyz = 1 and (xyz)" = 1 are "the same", 
care needs to be taken. If n is an even integer then the surface will have twice 
as many components as when n is odd; for example the point (1, 1,-1) is on 
x'^y'^z^ — 1 but not on xyz = 1. To show they are isometric would require the 
calculation of the first fundamental form, which is not well-defined for algebraic 
surfaces, i.e. surfaces defined implicitly by F(x,y,z) = c. We can calculate the 
Gauss curvature using the following expression [13| (here VF and H{F) are the 
gradient and Hessian of F respectively, and the norms are w.r.t. the ambient 
Euclidean space) 
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which we find to be independent of the value of n, but having the same Gauss 
curvature at identified points is a necessary but not sufficient condition for 
isometry. Instead, we can generate the geodesic equations themselves on the 
algebraic surfaces in question, and prove the following theorem. 

Theorem 2. The geodesic flow on the algebraic surface a;"?/"z" — 1 with n G M 
is not Louiville integrable with meromorphic first integrals. 

Proof. The geodesic equations on the algebraic surface F{r) ~ c where r = 
(x, y, z) are given by [6| 

{H{F)r).r 



Taking F — x^y'^z^ we find the geodesic equations are independent of n, i.e. 
the geodesic equations are the same for all values of n (we need to make use of 
the fact that F = VF.r — 0). We have shown the geodesic equations are not 
Louiville integrable when n = 1 in the previous Section, and therefore they are 
not integrable for n € M. D 

4-2. Surfaces of the form x"y"'z = 1 

It might be hoped that we could generalize the surface considered in Section 
3 to Monge patches of the form 

n e N. (8) 



Unfortunately the techniques described in this paper do not allow for a uniform 
treatment, for the following reason. 

Using the methods of Section 2 the NVE of the planar geodesic of ([5]) is 

„ 2n2(2n + l)(4n3 - lOn^ - y'^'^+^{An + 5)) 

^ + 2/2(4^2+ y4n+2)2 ^ " 0- 

As before, there are regular singular points at 0, oo and the An + 2 roots of 
4n^ + y4"+2 = which are distinct and symmetrically distributed along a circle 
centred on the origin of the complex plane. The /3 coefficients are 

^ f (l + 2n)(2n-5) 5 5_\ 

Pj-<^ 4 '16'---'16i' '^°°~ ' 

where ^ appears An + 2 times. We note that 



Vl + 4/3o = 2^Jn?^2n~l ^ Q V n G N. 

To show this we note that i/l + 4/3o e Q ^ n^ - 2n - 1 = m^ for some m e N. 
Since n>m^m = n — rj for r; S N which leads to a quadratic in rj whose 
roots are (2?! + 5,-5) ^ N. Therefore we can rule out cases / and /// of the 
Kovacic algorithm as in Section 3. 

However, in analysing case //, the Ej sets as in ([2]) are 

E„ ^ {2}, -Bi...(4„+2) = {2, 5,-1}, -Boo = {0, 2, 4} 

and as such there will be a combination leading to d = ^(4 — (2 — (An + 2))) = 
2n + 2 and all values below. Thus we can at best look at individual values of n, 
for example: 

Theorem 3. The geodesic flow on the surface x^y^z = 1 in not Louiville inte- 
grable with meromorphic first integrals. 



Proof. There are 615 combinations of the indicial exponents leading to d = 0; 
55 leading to d = 1, 2, 3 and 1 leading to d = 4, 5, 6. Each of these need to be 
check as described in Section 3. Again, the procedure is tedious rather than 
difficult. As there are no combinations for which the necessary P exists, we can 
rule out case // of Kovacic's algorithm. Thus the identity component of the 
differential Galois group of the normal variational equation is not abelian, and 
therefore the geodesic flow on the surface z = l/(x^ — y^)^ is not integrable. D 

As the number of cases which need to be checked increases rapidly for in- 
creasing n, the methods described in this paper are not appropriate for testing 
the integrability of surfaces of the form x'^y'^z = 1. Having said that, since the 
n — 1 and n — 2 cases are not integrable, we would conjecture that all surfaces 
of this form with n G N are also non-integrable. 

5. Conclusions 

Using Morales-Ramis theory and Kovacic's algorithm we are able to rigor- 
ously prove the (meromorphic, Louiville) non-integrability of the geodesic flow 
on certain algebraic surfaces. This approach is very geometrical in flavour, as 
opposed to the topological approach followed by Kozlov [Tj]; nonetheless it is 
robust enough to deal with free parameters and perturbations as another paper 
by the author has shown (l5i|. The analysis was facilitated by two features of 
the surfaces considered: Monge patches allow a simple intrinsic coordinate sys- 
tem/parameterization to be defined, and a plane of symmetry leads to a planar 
geodesic along which the variational equations decouple easily. It would be of 
interest to consider other surfaces where these properties do not hold. 
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